Currents in a many-particle parabolic quantum dot under a strong magnetic field 
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Currents in a few-electron parabolic quantum dot placed into a perpendicular magnetic field are 
considered. We show that traditional ways of investigating the Wigner crystallization by studying 
the charge density correlation function can be supplemented by the examination of the density- 
current correlator. However, care must be exercised when constructing the correct projection of the 
multi-dimensional wave function space. The interplay between the magnetic field and Euler-liquid- 
like behavior of the electron liquid gives rise to persistent and local currents in quantum dots. We 
demonstrate these phenomena by collating a quasi-classical theory valid in high magnetic fields and 
an exact numerical solution of the many-body problem. 

PACS numbers: 73.21. La, 71.10.-w, 75.75.+a 



I. INTRODUCTION 

The problem of the electronic structur e - 1 ' 2 ' 3 of few- 
electron quantum dots^ still remains in the center of 
attention of solid state research. Although recently 
substantial achievements have been made by relying 
on the exact numerical solution^Sil of the complicated 
quantum-mechanical problem, simplified approaches^ 
and simple analytical models are also of great interest 
and even become increasingly more popular AiSiiiiiLiii^ 
Besides their relative simplicity, these models are attrac- 
tive due to the provided physical insight and transparent 
visualization of the relevant phenomena. 

The successful development of the approximation 
based on the picture of a rotating electron molecul o 12 ' 13 
encouraged us to look at the rotation in quantum dots 
more closely. Namely, we investigated the distribution of 
currents in the simplest quantum dot — the rotating elec- 
tron ring formed in a few-electron system in a parabolic 
confinement and a perpendicular magnetic field. 

Bedanov and Peetersi^ predicted the crystallization of 
a system of classical two-dimensional (2D) particles in 
an external parabolic confinement into a set of concen- 
tric rings. In quantum dots containing up to five elec- 
trons only one ring is formed. An accurate quantum- 
mechanical solution based e. g. on the exact numerical 
diagonalization of the many-electron Hamiltonian pro- 
vides a circularly symmetric distribution of the charge 
density. The formation of the Wigner crystal can be seen 
in the density-density correlation functioni£*i£ obtained 
by considering the conditional probability to find an elec- 
tron at a point r given that there is an electron at another 
point ro. 

Recently, the analysis of various electron structures 
was further stimulated by employing a description based 
on mean-field approaches, such as the density-functional 
theoryj 3 " In this method, the electron-electron correlation 
is taken into account as an effective single-particle po- 
tential, and the correlation function is not easily acces- 
sible. Thus, the rotational symmetry of the dot has to 



be broken in some artificial way, and the electron crys- 
tallization at high magnetic fields is made visible directly 
in the electron density. Consideration of the currents^ 
also portrays the presence of a crystal — currents are 
circulating around the charge density lumps. 

In the present paper, we study the local currents ap- 
pearing in the vicinity of the electron lumps in a rotating 
Wigner molecule from the point of view of the density- 
current correlation function. This function describes the 
distribution of conditional currents in the system given 
that there is an electron at a certain point and indeed 
shows the electron crystallization. However, a certain 
care must be exercised when constructing 2D projections 
of the multi-dimensional space of many-body currents. 

The magnetic field tries to rotate the electron system 
as a rigid body, i. e. with a constant vorticity. As the 
Schrodinger equation has no dissipation the "quantum 
electron liquid" is akin to the Euler liquid in hydrody- 
namics and tries to rotate itself with a fixed angular 
momentum, that is, with vorticity equal to zero every- 
where except at the origin. This contradiction leads to 
two interesting phenomena: local currents and persistent 
currents. Both of them are caused and controlled by 
the electron-electron interaction. We demonstrate these 
phenomena by constructing a simple quasi-classical the- 
ory employing a rotating frame in the limiting case of 
high magnetic fields, and illustrate the obtained results 
by comparing them to the exact numerical solution of the 
considered problem. 

The outline of the paper is as follows. Section II de- 
scribes the model used in our calculations. The next 
two Sections are devoted to the analytical solution: in 
Sect. Ill we derive the Hamiltonian in the rotating frame, 
and Sect. IV describes its quantization. The results con- 
cerning the persistent and local currents are presented 
in Sects. V and VI, respectively. We conclude with a 
summarizing Section VII. 
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II. MODEL 

We consider a 2D parabolic quantum dot containing 
N = 2 — 5 electrons placed into a perpendicular mag- 
netic field B. This system is known to form a single 
electron ring. The magnetic field is described in terms of 
its symmetric-gauge vector potential A = |[B x r], then 
using the standard procedural of switching to dimension- 
less variables we obtain the following Hamiltonian 
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Here, the energy is measured in units hwo with ojq being 
the characteristic frequency of the confinement potential. 
The coordinates are measured in units Iq = y/h/m*u>o 
(the oscillator lengths), and the magnetic field in units 
$o/tt'o (^o = nhc/e). The dimensionless Coulomb cou- 
pling constant A = lo/a* B is expressed as the ratio of the 
parabolic confinement length Iq to the effective Bohr ra- 
dius a* B — eh 2 /m*e 2 . One way to find the eigenstates of 
the above Hamiltonian is to perform the "exact" numer- 
ical diagonalizations in the truncated Hilbert space of all 
possible electron configurations as described in e. g. Ref. 
5. Such calculations are feasible for few-electron quan- 
tum dots. Concentrating on the orbital effects we neglect 
the Zeeman energy but do take the degeneracy due to 
the spin states into account when constructing the ba- 
sis of many-body states. This approach is supplemented 
by an approximate analytical scheme valid at high mag- 
netic fields and employing a quasi-classical expansion in 
powers of B~ 1 whose development is the subject of the 
following Sections. 



III. ROTATING FRAME 

Since we are looking for the ground state wave func- 
tion, the strong magnetic field term in the parentheses of 
Eq. JP must be compensated by a large angular momen- 
tum M that appears in the first term of the same paren- 
theses when the gradient operator acts on the rotating 
electron ring wave function. The easiest way to take this 
cancelation of terms into account is to use the Eckardt 
framed which is the main instrument employed in the 
analysis of the rotation-vibration spectra of molecules 
and was recently adapted to quantum dotsii& 

The idea of the rotating Eckardt frame is to decouple 
the two degrees of freedom related to the center-of-mass 
motion as well as one more degree of freedom of the rota- 
tional motion of the system as a whole around its center 
of mass. The angular velocity of the rotating frame is 
chosen so that the system in the Eckardt frame has zero 
total angular momentum. As we are going to apply the 
expansion in B~ x powers and are interested only in the 
lowest-order harmonic approximation for the vibrations 
of the Wigner crystal, the above technique can be sim- 




FIG. 1: Local coordinates for three electrons on a ring. 



plificd using a frame rotating about the center of the 
quantum dot rather than the center of mass. 

In order to derive the Hamiltonian written in the ro- 
tating frame we follow the suggestion of Ref. |Ta and start 
from the classical Lagrangian 



L — Lmae — V. 
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where L mag is the sum of free-electron Lagrangians in 
a magnetic field, and V stands for the confinement and 
interaction potentials. 

In order to transform the above expressions into the 
rotating frame we follow Ref. |2(] and introduce the lo- 
cal coordinates (as shown in Fig. ^ for the case of three 
electrons). We denote the classical equilibrium radius of 
the ring by a. The azimuthal angles a n = ot(n — 1) with 
a = 2tt/N and n = 1, • • • , N indicate the equidistant lo- 
cations of electrons on the ring. We choose the local axes 
x n to be parallel to the equilibrium location vectors a„ of 
the respective electrons, and y„-axes are directed along 
the ring in the positive (counter-clockwise) direction. We 
assume that the whole frame shown in Fig. ^ is rotating 
in the positive direction with a constant angular velocity 
X- Thus, the position of the vector a„ is given by the 
angle 



a n (t) = a n + x{t), X{t) = Xt- 



(3) 



We represent the electron coordinates r = {x, y} by a 
complex number z = x + iy, and introduce the above- 
described local coordinates by means of the transforma- 
tion 

z n — > exp(i[a„ + x])(a + z n ). (4) 
This enables us to present the magnetic part of the La- 
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grangian (|2bll as 
1 N 

i mag = - + - fi2 /4) \a + z n \ 2 + 2ulm(z n z*)} , 

n=l 

(5) 

with u> = x ~ -B/2. This function is constrained by the 
condition 



/v 



/v 



Im ^ (a + Z n )z* w a Im ^ i* = 0, 



(0) 



which expresses the equality of the total angular momen- 
tum of the electron system to zero. Integrating the above 
constraint we rewrite it as 



JY 



E Vn = 0, 



(7) 



which can serve as the definition of the rotating frame. 

In order to obtain the corresponding Hamiltonian we 
introduce the normal modes taking into account the fact 
that the electron ring is invariant with respect to rota- 
tions by a multiple of a. Thus, the above normal modes 
are just the Fourier transforms w k — u k + iv k 
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Note that the index k numbering the normal modes runs 
from to N — 1. Inserting these expressions into Eq. (JSJ) 
we obtain the following final Lagrangian 

imag - \ {tio + - B 2 /4){^Na + u ) 2 } 
1 

+ o E il^l 2 + ( u * - s2 /4)kfc| 2 + 2ulm(w k w* k )} .(9) 
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The main advantage of this expression is that the condi- 
tion J7J is taken into account automatically by excluding 
the Vq mode which is now replaced by the rotation angle 
X- Next, we introduce the canonical momenta 
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where 



N-l 



Z = J +]TK| 2 , lo = {VNa + u ) 2 (11) 



is the moment of inertia of the electron ring. In order 
to obtain the standard form of the magnetic part of the 
Hamiltonian 



N-l 

E 
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Mx + U u +y(U k u k + V k v k )-L mag (12) 



we have to solve Eqs. (110(1 for velocities. However, aiming 
to arrive at the Hamiltonian in the harmonic approxima- 
tion we are entitled to make some approximations in the 
solution. Namely, in Eqs. I|10cfl and (|10d|) we replace 
the frequency by its approximate value u> ~ A1/2q 7 solve 
them for velocities ii k and v k , and inserting the obtained 
values into Eq. I|10a|) we obtain the following approximate 
expression for the angular velocity 



B M 1 %-^. r tt . 
X ~ — + — -— > (V k u k -U k v k )- 
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(13) 



Inserting these expressions of velocities together with the 
Lagrangian © into Eq. 1|12[) we arrive at the final expres- 
sion for the magnetic part of the classical Hamiltonian 




Mu k 



.(14) 



The total Hamiltonian is obtained by adding the poten- 
tial (|2c|) consisting of two terms that are small in the 
B — > oo limit. 



IV. QUANTIZATION 

The quantization of the ring rotation is trivial because 
the rotation angle x is a cyclic variable. Thus, the orbital 
momentum M is a constant of motion, and quantizing we 
simply replace this momentum by an integer eigenvalue 
of the corresponding angular momentum operator. The 
rotational part of the wave function is given by ^ ro t = 
exp(iMx). 

Before proceeding to the quantization of the remaining 
vibrational modes we have to fix the radius of the ring a. 
It is convenient to define this parameter by minimizing 
the potential energy which consists of the first term of 
the magnetic Hamiltonian (|14fl and the weak potential 
(|2cf) . Introducing the deviation A = X — Tq and replac- 
ing the electron coordinates in the second term of the 
weak potential by their equilibrium values, we present 
the potential as 



V = V - 



1 
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where the equilibrium potential is 
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k =l 
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and the factor /at = \ J2n=i I sin(an/2)| _1 is the 
Coulomb energy per electron in a ring of unit radius. 
Minimization of the equilibrium potential defines the 
equilibrium radius of the electron ring. As the first term 
in Eq. I|f 6|l is large and the second one is small, it is easier 
to obtain the estimate of the radius in two steps. First, 
we equate the derivative of the first term of the potential 
((Ttj)) to zero and find 
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M = \M\ = -Bl = 7;NBal 



(17) 



This expression gives us a relation between the magnetic 
field strength and the angular momentum in the ground 
state, and confirms the already mentioned fact that the 
angular momentum of the system grows in absolute value 
with increasing magnetic field. Then, the radius of the 
ring ao itself is obtained by minimizing the second (small) 
contribution to the potential 
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FIG. 2: The radial distribution of the density (p) and az- 
imuthal current density for three electrons in a quantum 
dot with A = 4 for angular momentum M = —18. The for- 
mation of a ring structure is apparent. Note that the currents 
flow in opposite directions on the inner and outer edges of the 
ring. The peak of p and zero of j v are close to the classical 



2 a 2 $ N ~ ^' (^) electron ring radius ao ~ 1.322 marked by the dot. 



with the result 



ao = (A/jv/2) 1 / 3 . 



(19) 



Now the construction of the Hamiltonian in the har- 
monic approximation can be completed. Taking into ac- 
count that the fulfillment of the condition (|17|l zeroes the 
second term of the expansion l|15|) , and using the approx- 
imation A ss 2ao\/N uq in the third term of Eq. i|15fl . we 
obtain the vibrational part of the Hamiltonian: 



B 2 ul) 





(20) 



The first term of this Hamiltonian (corresponding to 
the breathing mode) is just the Hamiltonian of a har- 
monic oscillator, and its ground-state eigenfunction is 
exp(— Buq/2). The other part of the Hamiltonian rep- 
resents a collection of nonintcracting two-component 
modes whose Hamiltonian resembles the Hamiltonian 
of a 2D electron in a perpendicular magnetic field. 
Thus, the corresponding ground-state eigenfunction is 

exp^EtYKlW 

Collecting all parts of the wave function together and 
performing the inverse Fourier transformation we obtain 
the electron ring ground state wave function 



K = 



e iM Xe -BK/i^ 

, 2 {J2n=l x n) (En=l 



(21a) 
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V. PERSISTENT CURRENTS 

Let us proceed to the examination of the current flow 
in quantum dots. In Fig. [21 the distributions of the 
charge and current density obtained from the exact di- 
agonalizations are plotted. Since these distributions are 
circularly symmetric, only the radial dependences need 
to be shown, and moreover, the current possesses only 
the azimuthal component j v . Fig. [3 is obtained for a 
three-electron quantum dot with the interaction constant 
A = 4. The angular momentum is M — — 18 and the di- 
mensionless magnetic field is set to B = 6.78, i. e. the 
midpoint of the magnetic field range where the consid- 
ered state is the ground state. 

We see that at these parameter values the electron ring 
is already well defined; the electron density at the origin 
drops almost to zero and peaks close to the classical value 
of the radius ao = (4/v / 3) 1 ^ 3 ~ 1.322. It is interesting to 
observe that the currents flow in the opposite direction 
on the inner and the outer circumference of the ring and 
vanish on the classical radius. 

Also, it can be noted already from Fig.Qlthat the radial 
dependence of the current density is nearly antisymmet- 
ric with respect to the classical radius, and thus, the two 
currents flowing in the opposite directions nearly cancel 
each other. That is, there is almost no global current 
running along the electron ring. In order to investigate 
this phenomenon in more detail we numerically calculate 
the magnetic field dependence of the net current crossing 
the dot radius 



/ = 



j v (r)dr 



(22) 



in the ground state. The results are shown in Fig. EI We 
consider angular momenta up to \M\ = 20, and use the 
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FIG. 3: The net current crossing the radius of quantum dots 
containing two, three and four electrons calculated by exact 
diagonalization. The angular momenta up to \M\ = 20 are 
included. We observe a saw-tooth-like behavior reminiscent 
of persistent currents in quantum rings. 



value A = 4 for two- and three-electron dots while in the 
case of four electrons in the dot we set A = 2. 

We see that these dependences display a saw-tooth-like 
behavior. The net current I increases continuously with 
the magnetic field in each ground state of a given angular 
momentum and drops abruptly at each angular momen- 
tum transition reminding persistent currents in a ring. 21 



These currents are responsible for the dot magnetization 
which demonstrates a similar behavior.— Note that in 
two- and four-electron quantum dots where the ground 
state at low magnetic field has the angular momentum 5 
AI = there is no net current at zero magnetic field. In 
contrast, the B — ground state angular momentum in 
three-electron quantum dotsi is \M\ = 1 and we observe 
a large current at low fields. 

For the sake of reference, we denote the absolute mag- 
nitude of the current jumps by AI and the average value 
of the current at the jump by /. In general, both AI and 
/ decrease monotonically towards zero with increasing 
magnetic field with a few notable exceptions occurring at 
very low values of the magnetic field. That is, the current 
oscillations tend to decrease in magnitude with increas- 
ing B and the oscillatory saw-tooth pattern becomes cen- 
tered around 1 = 0. The case of two electrons is rather 
regular and thus easy to analyse. Here, I reaches values 
very close to zero already at angular momenta around 
\M\ ps 10 while the behavior of three- and four-electron 
dots is more complicated. There are conspicuous irreg- 
ularities associated with more stable states at the magic 
values of angular momentum, and moreover, the oscilla- 
tions do not center around / = up to higher magnetic 
fields. For three-electron dots I approaches zero only at 
highest values \M\ « 20 while in four-electron dot the 
current oscillations do not center around zero in the con- 
sidered parameter range at all. 

We find that the dependence of I on the magnetic field 
is approximately exponential which indicates its essen- 
tially quantum-mechanical nature. Therefore, this as- 
pect of current oscillations can not be captured by our 
simple quasi-classical model whose prediction is I = 0. 
But in contrast, the behavior of AI can be analysed and 
understood from a classical point of view. 

For this purpose, it is sufficient to approximate the 
angle \ as 



1 



JV 



71=1 



(23) 



and the density of the electron current along the ring can 
be calculated as N times the average of the one-electron 
velocity operator. Its sole azimuthal component reads 



and the total current can be estimated as 
' M 



(24) 



I = 



dr 
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TjBNr^J p(r) j J d 2 rp(r) 
drp(r) / / d 2 rp(r) 
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BNa 



(25) 



Thus, the current is proportional to the largest factor in 
the potential (|16|l term which was assumed to be zero 
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when denning the approximate ring radius ao- It means 
that for the calculation of the persistent current the ra- 
dius of the ring has to be defined with a greater precision. 
So, let us equate the derivative of Eq. to zero 

d 1 ( 9 B 2 N 2 a i \ ( NXf N \ 

Using this expression together with Eq. (|17|) we present 
the current as 

7=^fM-WV 1 (a- »*!L) 

N f Xf N f NB\ 3/2 \ l2 ' ' 

2wB I 2 V 2 |M|y J' 

We observe that the magnetic field can not be compen- 
sated exactly by a discrete value of the orbital momen- 
tum. Thus, we define the critical values of the magnetic 
field Bm that correspond to the exact compensation as 

2\M\ 2\M\ , , , 9/q 

Then, introducing the magnetic field deviation B = 
Bm + AB we rewrite the expression for the current l(77|) 
as 

7 = -^{l-(^±^y /2 L^. (29) 
2ttB m \ V B M J J 4ttB^ v ' 

Finally, estimating the value of the current jump at the 
angular momentum transitions we substitute Ai?( max ) = 
2/Na 2 , into Eq. J22J and find 

A 3 
AI= ^k' AN= 2,(f N /2)^ (30) 

Thus, our theory predicts a power-law dependence of 
the current jumps on the magnetic field AI ~ _B~ 2 . In 
order to test this prediction against the results of the nu- 
merical calculations, in Fig.0]we plot AI as a function of 
B for quantum dots containing two and three electrons. 
The symbols depict the numerical results from our ex- 
act diagonalizations while the straight lines describe the 
quasi-classical limit. One sees, that at high magnetic 
fields the numerical results for two- and three-electron 
dots indeed cross into a linear regime in a full agreement 
with Eq. (jUJ. 

The distribution of currents in quantum dots shown in 
Fig. H and the appearance of saw-tooth-like oscillations 
can be explained on the basis of the following simple con- 
siderations. The physical current consists of two compo- 
nents, the so-called paramagnetic current proportional 
to the gradient of the phase of the wave function and 
the contribution due to the vector-potential. The former 
component resembles the irrotational flow of the Euler 
liquid known in hydrodynamics. This current behaves 




ln(B) 

FIG. 4: Double-logarithm plot of the current jumps A7 
versus the dimensionless magnetic field B. The symbols 
(full lines) denote the numerical exact diagonalization (quasi- 
classical) results. 



as j v oc m/r where m is the electron orbital momen- 
tum, and therefore, its vorticity (curl) is zero everywhere 
except at the origin. On the other hand, the vector- 
potential part introduces a rigid-body-like rotation with 
jip oc r, and the vorticity of this component is the same at 
every point. As any physical system tries to minimize its 
energy and tends to have as small currents as possible, 
the two components are required to cancel each other. 
However, due to the very different r-dependences such 
cancelation is possible only at the classical radius of the 
quantum dot. 

The net integrated current crossing the radius of the 
quantum dot is also, in general, non-zero. According to 
Eq. I|17l) , an exact cancelation of the net current can be 
realized only for a specific value (close to the classical 
value) of the electron ring radius. The Coulomb repul- 
sion between the electrons precludes the adjustment of 
the radius thereby causing the appearance of the uncom- 
pensated persistent current. This many-electron effect 
is essentially different from the single electron models^ 
where persistent currents appear due to fixation of the 
electron ring radius by the confinement potential or sys- 
tem boundaries. 



VI. LOCAL CURRENTS 

In order to obtain a better understanding of the com- 
peting currents flowing in the opposite directions on the 
inner and outer edges of the electron ring, we will now 
consider the density-current correlation functions which 
could also be called conditional currents. In close similar- 
ity to the more familiar density-density correlators (con- 
ditional densities), these functions are obtained by pin- 
ning one of the electrons at a certain point and inspecting 
the distribution of currents created by the other elec- 
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trons in the dot. In our calculations, we place the pinned 
electron at the distance equal to the classical radius do 
from the dot center. The remaining electrons will tend to 
localize close to their crystallization points, distributed 
cquidistantly along the circumference of the ring. In the 
quasi-classical wave function l|21|) the electrons are distin- 
guishable and the motion of each of them is restricted to 
the vicinity of its own crystallization point. Therefore, 
when evaluating the quasi-classical correlation function 
close to these points we may dispose with the summa- 
tion over the electrons and take into account only the 
nearest one. 

When considering the correlation functions, the ro- 
tating frame coordinates (x n ,y n ) are not the most con- 
venient ones due to the presence of the constraint JJJ). 
Thus, having a pinned electron does not imply fixed val- 
ues of its coordinates. In order to get rid of this nuisance 
we switch to a new set of coordinates rj n ) which mea- 
sure, respectively, the radial and angular deviations of 
the electrons from their crystallization positions in the 
laboratory frame. In the considered complex representa- 
tion {2J both coordinate sets are related as 



exp(i[a„ + x])(a + x n + iy„ 
= exp(ia„)(a + £„ + ir? n ). 



(31) 



Separating the real and imaginary parts we are immedi- 
ately led to the following two equations of the coordinate 
transformation 



(£„ + a) cos x + Tin sin x, (32a) 
?7„cosx- (£« +a)sinx- (32b) 



fjn 



Performing the summation over all electrons in Eq. Il32bl) 
and fulfilling the constraint Q we define the rotation 
angle \ 



tanx 



Na + Y,n=l £™ 



(33) 



Now, inserting this expression into Eqs. 1-S-?I) we solve 
them for (x n ,y n ). Further inserting the solutions into 
Eq. (|21|l we obtain the quasi-classical wave function in the 
laboratory frame. Taking into account the fact that the 
quasi-classical wave function was derived in the quadratic 
approximation and correspondingly retaining only the 
necessary terms 



T N n 



1 - 



Na 



1 N 



(34a) 
(34b) 
(34c) 



we obtain the wave function in the laboratory frame as 
* = e iMx e -BK/i (35) 



Here, the function K is defined by the same expression 
as Eq. (|21b|) with the variables (x n ,y n ) replaced by the 
laboratory- frame coordinates (£ n , r) n ), and the angle x is 
now given by Eq. (|34a|l . 

We begin by investigating the shape of the electron 
density lumps formed in the vicinity of their crystalliza- 
tion points as these results will turn out handy when 
discussing the currents. Straightforward integration of 
the iV-particle density 



Pn 



1*1 



-BK/2 



(36) 

over the coordinates of N — 2 electrons with the N th 
electron fixed (i. e., £/v = T)n = 0) and the first electron 
coordinates renamed into £, 77, gives 



exp 



d£ 2 di] 2 ■ ■ ■ d^N-idriN-i Pn 
B (2N- 1 



4 V N- 1 



>r 



(37) 



That is, the electron lump has the form of a Gaussian 
elongated in the azimuthal direction. The contour lines 
connecting the constant-density points are ellipses with 
the ratio of semiaxes 



Pn = 



2N- 



N - 1 



i = V3, 



for N = 2,3,4. 



(38) 

We calculate the current-density correlation function 
in two steps. First, we consider the total current in the 
TV-particle space and write out its component due to the 
first electron 



ji,jv = {MVix + A}/9jv 



M 

.zip) , 
' Jl,jV ^ 



(iVa )2 

M En=l 



B 



B 



(Na y 



+ ir( a o + Ci) 



(39) 



>Pn 



.5(0 
Jl.JV 



Here, the symbols and denote the unit coordinate 
vectors. The term independent of £ and 77 



3l > N ~ \Na Q 



Ba 



e n p N 



(40) 



leads to the persistent current which was already dis- 
cussed above. Due to the integer-valued quantization of 
M this current vanishes only at the special values of the 
magnetic field strength Bm given by Eq. I|28|l . Using 
Eq. H17|) the term linear in £ and rj describing the local 
currents can be rewritten as 



:(0 
Jl.jV 



B 



1 



N 
n=l 



1 N 



(41) 



>PN 



x Vi]pjv, 
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where e 2 is the unit vector along the direction of the 
magnetic field. Note that this relation has the same form 
as that previously obtained for single-electron densities 
and currents^ 

Now we are ready to perform the second step and cal- 
culate the correlator of the local currents by integrating 
the obtained expression l|41|l over the coordinates of N— 2 
electrons and pinning the ./V th electron as it was done in 
obtaining Eq. J3JJ|. This procedure amounts to the re- 
placement of the ./V-particle density pat by the density 
correlator l|37|l because the integration does not involve 
the coordinates of the remaining first electron. Thus, 
dropping the indices and denoting the first electron co- 
ordinates by (£, rj) we arrive at the final expression for 
the current-density correlation function 



x V]p(£,r/). 



(42) 



The obtained simple expression leads to important 
consequences. The current lines are perpendicular to the 
gradient of the density, and 



divj(£,??) 



1 



V[e 2 xV]/)K,i))=0. 



(43) 



That is, in the considered quasi-classical approximation 
the local currents circulating around the localized elec- 
trons are conserved. Therefore, they are physically well 
defined even though there is no general conservation the- 
orem for the conditional currents (see Appendix |A"|) . 

This circulatory motion of electrons is similar to the 
cyclotron rotation present in any single-electron model 
(see, e. g., Ref. 0). However, due to the electron corre- 
lation and according to Eq. I|42[l the electrons rotate not 
along Larmor circles but along elliptic density contour 
lines. The standard quantity describing such rotational 
motion is the vorticity (curl) of the current field. In our 
case it reads 

curlj(£,»?) = --[V x [e, x V]]p(£,?/) = -iv 2 p(£^) 



= j{Pn 



+ i-B( P %e+v 2 )}p^,v)- 



(44) 



It is worth pointing out that the vorticity is positive (i. e., 
it has the same sign as the vorticity of the current com- 
ponent due to the vector potential) at the electron crys- 
tallization points, and becomes negative at a certain dis- 
tance away from them. The constant-vorticity contours 
are more elongated in the azimuthal direction than the 
density contour lines, in particular, the contour corre- 
sponding to the zero vorticity has an elliptic shape with 
the ratio of semiaxes equal to p 2 n rather than p n . 

We illustrate the above conclusions with numerical re- 
sults. The divergence of the density-current correlation 
function calculated for a two-electron dot is plotted in 
Fig. We choose the value A = 4 so that the classical 
quantum dot radius (indicated by a dotted line in the 



>- 




0.0 0.5 1.0 1.5 2.0 



X 





0.0 0.5 1.0 1.5 2.0 



FIG. 5: The divergence of the density-current correlator in a 
two-electron quantum dot. The angular momentum is M = 
— 18. Panel (a) corresponds to the the lower end (i. e., B = 
17.75) of the magnetic field interval where this state is the 
ground state, panel (b) is plotted for the medium value (B = 
18.25), and panel (c) is obtained at the higher end (B — 18.75) 
of this interval. Uncompensated global currents are visible. 
Dark (light) areas correspond to generation (extinction) of 
the current. The dotted curve denotes the classical radius of 
the quantum dot. 
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plots) is do = 1, and place the pinned electron at (—1, 0). 
Dark (light) areas correspond to the positive (negative) 
divergence, i. e. generation (extinction) of the current. 
The three panels of Fig. [5] show the same ground state 
with \M\ = 18, however, at slightly different magnetic 
fields. The plot in panel (a) is obtained for the lowest 
possible magnetic field at which this state is the ground 
state (B = 17.75), panel (b) corresponds to the mid- 
dle of this interval of magnetic fields {B = 18.25), and 
panel (c) is obtained at the higher limit of this range 
(B = 18.75). In all three plots non-conserved currents 
are visible. Panels (a) and (c) show a dipole-like struc- 
ture corresponding to the azimuthal persistent current 
1)41) [I. At low magnetic fields [panel (a)] there is a net 
current running in the clockwise direction and at high 
magnetic fields we observe a net counter-clockwise cur- 
rent, in accordance with Eq. (|40|l and the saw-tooth-like 
behavior of Fig.[3f a). The data in panel (b) is obtained at 
the middle of the allowed interval of the magnetic fields. 
Here, B = Bm and there is no global current. However, 
the divergence assumes a quadrupole-like checkerboard 
pattern still indicating the presence of a certain current 
non-conservation. This is a small higher-order effect not 
captured by the above quasi-classical treatment. 

Fig. EI a) shows the curl (vorticity) of a two-electron 
quantum dot. We find that the absolute values of the 
curl are 3 orders of magnitude larger than the diver- 
gence, so the current non-conservation is small indeed 
and can be safely neglected. Here, only one plot corre- 
sponding to the middle of the allowed energy range [as 
in Fig. OJb)] is shown since these plots are rather insen- 
sitive to small variations of the magnetic field. The first 
electron is pinned at the point (—1,0) on the classical 
radius (the dotted line), and the second electron is crys- 
tallized in the vicinity of the opposite point (1,0) and 
performs a cyclotron-resonance-like motion. Dark (light) 
areas correspond to the positive (negative) vorticity, and 
the full black line separates the areas corresponding to 
different signs of the curl. According to the theory pre- 
dictions, this line is an elongated ellipsis and is marked 
by a white line. We see that the numerical results indi- 
cate a certain spreading and deformation of this area as 
the quasi-classical regime is not truly reached. 

For three and four electrons we also obtain a dipole-like 
structure in the divergence plots since at the considered 
magnetic fields and angular momenta the global rotation 
of the electron molecule is not fully stopped, i. e. the os- 
cillations of the current [see Figs.|3Ib) and|3c)] are not 
centered around zero value. In these cases, the maximum 
absolute values of the divergence are around 25 times 
smaller than the corresponding magnitudes of the curl. 
The plots of the curl are shown in Figs.^b) and[|J{c) for 
three and four electrons, respectively. We again consider 
the states with the angular momentum M — —18, and set 
the magnetic field strength to the midpoint value of the 
field range corresponding to the considered ground state. 
One electron is pinned on the negative part of the x-axis 
at the distance equal to the classical radius (marked by 




-1.0 - 



-1.0 -0.5 0.0 0.5 1.0 

X 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



X 



2.0 




-2.0 I I 

-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



X 

FIG. 6: The curl of the density-current correlation function 
for quantum dots containing two, three and four electrons. 
The ring of the classical radius is marked by the dotted line, 
and one electron is pinned at the intersection of this ring 
with the negative part of the x-axis (indicated by the dot). 
The ground state of M = — 18 is shown. Dark (light) areas 
correspond to positive (negative) vorticity, and the full black 
line denotes the separating zero-vorticity contour. The quasi- 
classical prediction for this contour is marked by the white 
line. 
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the dotted line) from the center. For three electrons with 
A = 4 this radius is a — 1.322, and for the four-electron 
dot with A = 2 we have ao — 1.242. The plots in Figs. 
|SJb) and|SJc) again show the remaining electrons local- 
izing close to their crystallization positions and perform- 
ing the cyclotron-like motion along trajectories elongated 
in the azimuthal direction. The numerically calculated 
zero-vorticity contours are expanded and distorted with 
respect to the ones predicted by the quasi-classical ap- 
proximation. In the case of four electrons in a dot the 
regions of positive vorticity even overlap. 



VII. CONCLUSION 

In conclusion, we developed a quasi-classical theory 
based on a transformation to a rotating frame that is 
capable of describing the currents in quantum dots at 
high magnetic fields. The results show that due to 
the competition between the paramagnetic and vector- 
potential components, there arise global persistent cur- 
rents running along the electron ring. At the same time, 
each electron may be visualized as performing individual 
cyclotron-like motion along elliptic trajectories elongated 
in the azimuthal direction. This motion is well described 
by the density-current correlation functions. We intro- 
duced these functions as a refinement of the suggestion 
to demonstrate the Wigner crystallization by calculat- 
ing ordinary currents in the broken-symmetry situation, 
and showed that in the quasi-classical limit they are con- 
served, and thus, physically well defined. 
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APPENDIX A: TWO ELECTRONS ON A ID 
RING 

In this section, we present a simple solvable model that, 
by means of a geometrical interpretation, illustrates the 
problems arising with the single-electron functions and 
correlators when discussing Wigner crystallization. 

We consider two electrons moving on a ID ring of ra- 
dius r = 1 in a perpendicular magnetic field whose be- 
havior is described by the following Hamiltonian 

(Al) 




FIG. 7: Many-electron density. 



Here, the positions of the electrons on the ring are given 
by the angles <fi : 2, the symbol 7 = eBr 2 /2ch stands 
for the dimensionless vector potential, and the electron- 
electron interaction is modeled by the cosine function 
which ensures the Wigner crystallization of the electrons 
on the opposite ends of the diameter in the strong inter- 
action (A — > 00) limit. 
After the transformation 

$ = i(( / 9 1 +(^ 2 ), <f = (fi! - Lp 2 + IT (A2) 

the variables separate, and we rewrite the Hamiltonian 
(|A1|) in the A — > 00 limit as 

B --k{Si + »<)' (A3) 

whose ground-state eigenfunction is 

V(<Pl,<P2) = e iM (vi+f2)/2 e -( Vl - V2+7 r) 2 VV8. (A4) 

Here, the symbol M stands for the total angular momen- 
tum. In the ground state, its value equals to the integer 
closest to —2"f. The corresponding many-electron density 

p{i Pl ^ 2 )=e-^-^+^^, (A5) 

is shown in Fig. by the dark stripes in the 2D many- 
electron space. In this space, the current has two com- 
ponents related to both electrons. Its first component is 
given by the expectation value of the first electron veloc- 
ity operator (—id/dcpi + 7) 

j 1 = (M/2 + -y)p(<p 1 ,<p 2 ). (A6) 

The same result is obtained for the second component of 
the total current j'2. Thus, the current flows along the 
dark stripes as indicated by arrows in Fig. 



11 



Eqs. l|A5jl . ljA6|) and Fig. contain all available infor- 
mation about the system. However, usually it is too dif- 
ficult to consider many-electron spaces, and one resorts 
to simpler single-electron functions such as the density 
and the current obtained by integrating of Eqs. I|A5() and 
(|A6(I over the coordinates of the second electron 

= / d(p'p(ip, (p r ) = const, (A7a) 
Jo 

j(<p) = (M/2 + 7) / dcp'p(cp,<p') 
Jo 

oc M/2 + 7. (A7b) 

The current expression gives the persistent current flow- 
ing along the ring. From the geometrical point of view, 
the above expressions are projections of the 2D many- 
electron density onto the abscissa axis. Naturally, they 
lead to a homogeneous single-electron density which pro- 
vides no information about the Wigner crystal. Thus, 
one has to consider more sophisticated correlation func- 
tions. In our case, the density-density correlation func- 



tion is given by Eq. (|A5|) and the density-current cor- 
relation function by Eq. I|A6|I . In both cases, they 
are treated as functions of the first electron coordinate 
with the second electron coordinate fixed at the point 
if 2 = p' ■ Geometrically this actually means taking the 
cross-section of the plot shown in Fig.0along the dotted 
line ip' — const. The density-density correlation function 
is shown by a dashed line, and we see a lump at the point 
pi = ip = ip' + 7r which serves as a good indication of the 
formation of a Wigner crystal. 

The current-density correlation function is not so for- 
tunate because of the non-zero divergence 



-j-h(<P,<p') = -VA72(Af/2 + 7)(vJ-^+7r)p(^^) + 0, 

(A8) 

indicating the absence of the charge conservation. It is 
evident from Fig. [7| that the global current does not flow 
along the chosen cross-section, and the currents are es- 
caping into other dimensions. 
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